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Abstract 

A theory of double magnetophonon resonance (MPR) in quantum wells in 
magnetic field is developed. The magnetic field is assumed to be tilted at an 
angle 9 to the perpendicular to the plane of quantum well. The resonance is 
due to the resonant interaction of 2D conduction electrons with the longitu- 
dinal optic phonons. The electrons are assumed to be nondegenerate. The 
0-dependence of MPR maxima is investigated. The existence of a double reso- 
nance, i.e. two resonant peaks for each value of AA (the number of resonance) 
and the ^-dependence of the MPR maxima is explained by the screening in 
conjunction with the combined role of the phonon and electron damping and 
variation of the electron concentration in the well with magnetic field. 
PACS 73.50.Jt, 73.50.Mx, 63.20.Pw, 71.70.Di 

I. INTRODUCTION 

Magnetophonon resonance (MPR) in semiconductors is an internal resonance that is 
reached when the hmiting frequency of a longitudinal optic phonon equals the cyclotron 
frequency of an electron, Q, times some small integer. A/" (for instance, 0]). Since its 
theoretical prediction [0] and subsequent experimental observation in 3D semiconductor 
structures MPR, along with cyclotron resonance, has become one of the main instruments 
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of semiconducting compound spectroscopy. 

The advances in semiconductor nano-fabrication in recent years have made available 
materials of great crystalline perfection and purity. The electrical conduction and some 
other transport phenomena in such nanoscale structures has been a focus of numerous in- 
vestigations, both theoretical and experimental, with a number of important discoveries. In 
particular, the discovery of MPR in the quantum wells took place in the pioneering paper 
by Tsui, Englert, Cho and Gossard After this first publication a number of papers 
has appeared where various aspects of this physical phenomenon have been investigated. 
The most detailed experimental investigation has been done by Nicholas with co-workers 
(see the review paper P] and the references therein). It has been shown that the 2D MPR 
qualitatively differs from the same phenomenon in the 3D structures. As first steps in devel- 
oping theory of 2D MPR in the perpendicular magnetic field we can quote the theoretical 
papers They consider the MPR within the perturbation theory approximation with 

regard of coUisional broadening of the electron state. A theory in tilted magnetic field is 
developed in Ref. where transitions to a higher band of spatial quantization are taken 
into consideration. To our opinion, these theories do not provide detailed interpretation of 
the existing experimental data. Our purpose is to give an interpretation of experimental 
findings , such as the double resonance, by simultaneously taking into account the screen- 
ing of the phonon potential by 2D conduction electrons as well as the phonon and electron 



damping. For such a program, as we have indicated in [|T0| , the lowest approximation of the 
perturbation theory is not sufficient. As is shown in the present paper, the relative role of 
screening is determined not by the temperature (as has been suggested in Ref. |ll|) but by 
the interplay between the screening and the electron and phonon coUisional damping that 
is usually weakly temperature dependent. 

There are two main groups of the MPR experiments in quantum wells. The first group 
deals with the MPR in the perpendicular (to the plane of 2DEG) magnetic field. The main 
features of the findings in this case are (i) the fact that the resonance is determined by the 
transverse optic frequency ujt (rather than the longitudinal frequency uJi) and (m) a rather 



narrow interval of electron concentrations where the MPR is observable. The second group 
concerns with the experiments in a magnetic field tilted at an angle 6 to the perpendicular. 
Its characteristic feature may be called a double resonance. For small values of 9 the MPR 
is determined by uJt. Then, for somewhat larger values of 9 its amplitude sharply goes down 
within a narrow angular interval typically of the order of 10°. For even bigger values of 9 



there is another maximum, this time determined by uji [T^. These two types of resonance 
may be called the uot- and cj^-resonances as their positions are determined by the frequencies 
uit and uJi respectively. 

In our paper |]lOl we have given interpretation for the first group of experiments. In the 



present paper we offer interpretation of the second group. It is shown that the angular and 
concentration dependencies of the MPR amplitudes are deeply interrelated. 

The magnetic field B is assumed to be in the (y, z)-plane, the 2;-axis being perpendicular 
to the 2DEG, while the external electric field is oriented along the y-axis. The Pyy{= Pxx) 
component of the resistivity tensor will be calculated. This is the transport coefficient 
expressed through the 2D conductivity a^y (averaged over the width of the well) as 

Pyy ^xxl^^xy) ■ 

As usual, it is assumed that a1.y ^ o-^^o'yy. We consider the situation where the well is so 
narrow that only one electron band of spatial quantization is filled. Hence one can assume 
that the z-component of electron velocity vanishes. Then 

eB eB 
Px = eEx -\ VyCos9, Py = eEy VxCos9 (1.1) 

where px, Py and Vx Vy are respectively the components of the electron quasimomentum and 
velocity. These classical equations illustrate the physics describing a 2D circular motion of 
an electron with the angular frequency 

eB 

ncos9 = — cos^ (1.2) 
mc 

where m is the effective mass, so that p^^y = mVx^y It follows from Eq. ( |1 . Ij) that 



where e and n are the electron charge and concentration. It means that within one miniband 
approximation the magnetic field enters the problem only in the combination B cos 6. 



II. GENERAL EQUATIONS 

To calculate the x-component of the d.c. current it is convenient to consider the motion 
of a center of Landau oscillator. The conductivity axx averaged over the width of the well 
is given by (see Ref. JTB[ ) 

where S is the area of the 2DEG, T is the temperature, X is the operator of coordinate of 
the center of Landau oscillator in the Heisenberg representation. According to Eq. ( |1 . 1|) , in 
the Schrodinger representation 

ic d 

^ = n ^7r + ^ 

en cos u oy 

{cf. with Sec. 112). It commutes with the free electron Hamiltonian Ti in magnetic field 
B as well as with the operator of Coulomb electron-electron interaction. This is a conse- 
quence of the quasimomentum conservation in electron-electron collisions. Here (^X{0)X{t) 
is the ensemble averaged correlation function of velocities of the centers of Landau oscilla- 
tors. In the present and the following sections we will usually assume h = 1, ks = 1 and 
will restore these symbols in the resulting formulas. 
Now 

m = Y.j V^Hr,a)z[7^,X]^(r,a)rfV = ^-^^| i,\v,a)^i,{T,a)d\. (2.2) 

The summation is over the spin variable. Here ip is the operator of the electron wave function 
while U {t, r) is the operator of phonon field interacting with electrons. For the time being, 
we consider it as an external random field; later on we will average over all its realizations 



introducing the optic phonons. The expression for be presented in such a form (we 

remind that we calculate the conductivity averaged over the width of the well) 



a- 
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J dtj d^r J d^r' / n(0, r') ^f/(0, r')n{t, r)—U{t, r) \ (2.3) 



2TS\eBcosej 7-oo 7 J \ ' ' ' dy' ' ' ' ^ ' dy 
where n(t, r) the electron density operator. Representing the ensemble average in Eq. ( p.3| 
as a sum over the exact quantum states of the system (see Ref. |T^, Sec. 36) we get [jl6| 



axx = ^ r^^V f-f^fdz Idz'- 'y""^^^ 



2T\eBcoseJ J 2n J {271 f J J l-exp(-cu/T) 
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X [Dr{-u- q; z, z') - Da{-uj; q; z, z')] [uf {u; q; z', z) - lif {u- q; z' , z)] (2.4) 



(see the details of the derivation in Ref. [|T0[)- We made use of the quasimomentum conser- 
vation along the plane of the quantum well; q is a 2D wave vector parallel to the plane of 

/ON 

the well. N{ijo) is the Bose function. Iljj (cu; q; z', 2) is the exact 3D electron polarization 
operator. 

Now, Dr{uj\ q; z^ z') is the phonon propagator with regard of the direct Coulomb electron- 
electron (e-e) interaction [see below Eq. (|2.7|) ]. In the present paper we consider the mag- 
netophonon resonant contribution to axx- This means that the phonon contribution to the 
Green function is determined by the optic phonons. Further on we will assume that 
one can neglect the difference between the lattice properties within and outside the well. 
This assumption should not affect the qualitative results of the theory. [Eq. ( ^.41 ) permits 
to consider also a more general (nonhomogeneous) case]. Without regard of e-e interaction 
the phonon propagator has the form 



where k is the z-component of the wave vector while Ec is given by 

111 , , 

- = . (2.6) 

^00 ^0 

Here £0 and £00 are the dielectric susceptibilities for ^ and u ^ 00, respectively. We 
have included the Frolich electron-phonon interaction [|r^ into the definition of the zero-order 
phonon Green function. 



When calculating the exact phonon propagator it will be necessary to insert along with 
the phonon lines D^^\ the direct Coulomb interaction lines 

One should, however, observe the following important point. Both ends of the exact phonon 
propagator should be ordinary phonon lines D^^^ without Coulomb interaction lines. This is 
due to the fact that the operator X commutes with the electron-electron interaction operator 
[see Eq. ( p.2|) ] as the latter conserves the electron quasimomentum. 

Further we assume for the electrons a parabolic confining potential muj^z^ /2 with the 
following gauge for the vector potential A = {—BycosO + Bz sm6, 0,0). It is also assumed 
that 

huo > hQ, ksT (2.8) 

(where Q = eB/mc is the cyclotron frequency). In other words, we assume the confining 
potential to be strong. The diagonalization of quadratic Hamiltonian is a well-known pro- 
cedure (for instance, Ref. |]TB[). We wish, however, to emphasize that the actual form of the 



confining potential is not essential for the present theory provided that hui is much smaller 
than the distance to the bottom of the second miniband [see below — Eq. ( p.l2|) ]. 

The energy of confined electron in the magnetic field defined by the vector potential A 

is 

U = -mQ^ cos^ e{y - yof - -mQ^ sm{2e)z{y - yo) + -m(u;^ + Q'^ sin^ e)z^ (2.9) 
2 2 2 

where yo = —cpx/eBcos6 while px is the electron quasimomentum component that is con- 
served. We will see in Appendix that in the leading order in {Q/ujq)^ one can retain in 
Eq. (|2.9|) only the terms describing the electron motion in the magnetic field B cos 6 perpen- 



dicular to the plane of the well (c/. with Ref. |T9|)- This can be visualized in the following 



way. One can obviously neglect the magnetic field correction to the confinement potential, 
i.e. sin^ 9 as compared to Uq. This means that the characteristic values of z are of the 



order of / = ^Ujuj^m. The mixed term, i.e. the second term on the right-hand side of 
Eq. ( |2.9| ) for the typical values of z can be also discarded provided that fiuj^ is the biggest 
energy in our problem. Therefore the 2D polarization operator (calculated in detail in Ap- 
pendix) for the case of Boltzmann statistics and a small gas parameter we are interested in 
has the same structure as the polarization operator of Ref. [10| in the perpendicular field B 
with the replacement B Bcos6: 



Ur = -2ns exp 



{qas Y cotha 
2cos^ 



^ sinhATa ^ / q^o^l \ .3 10) 

^ " UVLcose + i5 1 2cos^sinha ' ■ ^' ' 



M=-oo 



Here Ij^ is the modified Bessel function, a = hVt cos 9 /2k-QT , a\ = ch/eB, Ug is the 2D 
electron concentration. The polarization operator n*-^-' of Eq. ( p.4| ) differs from Ur by the 
factor tp{z)'^ip{z')'^ due to the electron motion along the z-axis. Here iplz) is the wave function 
of the lowest level of transverse quantization. 

The e-e interaction can take place both via exchange of a phonon and as a direct inter- 
action described by V^'^\ Eq. (|2.7| ). The sum of two interactions is 

/ ,x 47re^ , , ujf — (uj ± iS)"^ . . 

Here uj^ =uf{l- £00 /^c) • 

Eq. ( |2.10D shows that the electron-electron interaction cannot be treated within the 
perturbation theory. Let the frequency u in Eq. ( p.lO|) be close to the frequency AfQ cos 6 so 



that only one term of the series is important. The higher orders of the perturbation theory 
(without regard of the electron damping Fg) give powers of an extra factor l/luj—N'Q cos + 
i6). Therefore, the e-e interaction must be a sum of chains of loop diagrams (see Ref. [0). 
Physically this means taking into account the screening of phonon polarization potential 
by the conduction electrons. Thus in 2D case in a resonance the screening can be very 
important. The reason as to why one does not need to take the screening into account in 3D 



case has been discussed in Ref. |T0[ . The only point demanding some attention is taking into 
account the spatial nonhomogeneity. However, the procedure is essentially facilitated by the 
fact that Il{zi, Z2) depends on il){ziYip{z2Y as factors. Then the index of the progression 
generated by the loops is proportional to 



dzidz2VR{q, 21 - Z2)ij\zi)ij\z2) (2.12) 



that can be presented in the form 



ds 



where ips is the Fourier component of ip{z). This resuh is vahd for a well of any form so far 
as the distance to the bottom of the second miniband remains much bigger than hui. It is 
only necessary to insert into Eq. (|2.12|) the appropriate wave function ip{z). 

For a quadratic confining potential we get a theory of 2D electrons with the interaction 
potential 

, , 47re^ f°° (ifcexp(-P/2/2) 
Further on we will also need the expression 

as well as the equation for the exact phonon Green function 

Vr[u}, q = — — -. 2.15 

1 + Vr{uj, q)nij(u;, q) 

Now we take into account the aforementioned point that both ends of the chain in 
Eq. ( p.4|) should be ordinary phonon lines (without the Coulomb interaction). We have 

DK(o;,q) =Di?)(.;,q) + Dg)(.;,q)— J- .Dt\u^,q) (2.16) 

where 

Dg)(u;,q) = j ^D^;^\u;,ci,k)eM-k'lV2) (2.17) 

and a purely 2D equation for the MPR 

^ 7 duj f d^q q^N{iu) 



a 



XX 



AT \B cose J J 2tt J (27r)2 1 -exp(-cj/T) 

— 00 

X [DR{-uj;q) - DA{-uj;q)][UR{uj;ci) - UA{uj;q)]. (2.18) 
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Here we made use of the fact that for huo ^ HQ the integrand in Eq. (|2.18| ) is symmetric 
in and qy. 

We will not insert directly Eq. ( |2.16| ) into Eq. ( |2.18| ) as it seems to have poles at u = ui 



that in fact disappear after integration and some algebra. It is convenient instead to present 
Eq. (^TTBp in the form 



Z..(.,q)^P.(.,q)-V(-.4^^-j^ 



(2.19) 



One can see that neither of these terms has a pole uj = uji [see Eq. ( |2.13| ) in combination 
with Eq. (PID I. This is a manifestation of the influence of screening. It means that the 



screening may play a certain role even for relatively small electron concentrations. 

The MPR is, as we will see, determined by the last term while the contribution of all 
the rest terms in Eq. ( p.l8 ) vanishes (provided that one neglects the electron and phonon 



damping). As a result, we have for the A/'th resonance of 

2nsc'^h'^ N(ut) smh(hiUt/2kBT) ^ ^.^ , , . ,^ 

= —n^2 — ^ 1 f \ /. 1.^ JmS[i e(mcose)] 2.20 

Sook^TB^ cos^ 6 1 — exp{—huJt/kBT) 

where 

= r i,iv'-\,)I^ { ^4^) exp . ,2.21) 

JO \ 2 cos 6 smh a J \ 2 cos 6 J 

As is indicated in Sec. |, it is natural that in the lowest approximation in Vt^ /uJq only the 
combination Bcos6 enters the equations describing the 2D motion of an electron in the 
quantum well. The integrals in the rest terms in Eq. (|2.19|) are either real or proportional 



to 6^^6{e~^). This can be easily checked if one takes into consideration that S{e^^) comes 
from the pole [11^^ + V/j(q)]^^ and the terms in question in Eq. (|2.19| ) have the factor e~^. 

Further we will be interested in the case ql <^ 1 when the effective e-e interaction does 
not depend on the form of potential and is equal to 

V(c)(q) = ^. (2.22) 

Then one can present as 



where Qjv'_i/2(2;) is the associated Legendre function of the second kind (we remind that 
Qa/'-i/2 is an imaginary function of a real argument). We will consider the case a -C 1. 
Then the characteristic values of q are of the order of qr = \/2mk-QT /h and one can present 
Jj^ in the following form 



As e'^^{uj) [Eq. has a zero at u; = Ut, a^x exhibits in this approximation an 

infinitely narrow magnetophonon resonance at 

Afn cos 9 = ujt. (2.25) 

Physically this is due to the fact that the e-e interaction without regard of the damping is 
very strong in the resonance. 



III. ANGULAR DEPENDENCE OF MPR MAXIMA 

In the present section we will investigate dependence of the positions of the MPR maxima 
on the angle 6. As we have indicated, in the limit of vanishing phonon and electron damping 
(r and Fe, respectively) the screening in the resonance is very strong. If we take into 
account that the dampings are finite one can calculate the critical concentration where 
the screening ceases to play a role. As the interaction depends also on 6 for each value of 
Us one can indicate the corresponding critical value(s) of 6. 

We will start with taking into account the phonon damping. Finite optic phonon damping 
is due to the decay of an optic phonon into two acoustic ones (see [|10|)- Technically it can 
be taken into account by replacement u u ± iT in the retarded and advanced phonon 
Green functions, respectively. One can easily see that in such a case the MPR acquires a 
finite width which one can take into account by the following replacement in Eq. ( ^.2(j| ) 
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Here 



Sie-')^-lmeR. (3.1) 

TT 



. (3^2) 



AfQcosO - Jlu! + T^) +r2 



TT 27r(jJt 

In what follows we will assume that 

Te, r < CJ; - U;t < UJt. (3.3) 

One needs the first inequality to be able to discriminate between frequencies ui and ut. The 
second inequality is fulfilled for such systems as GaAs/GaAlAs and facilitates the calcula- 
tions. 

Now we will discuss the role of the electron damping. Good examples of importance 
of the damping for the MPR are given in Refs. P,|20[]. For us it is important as it may 
both destroy the resonance and determine the angular interval for its existence we are 
looking for. We assume that Fe ^ QcosO. The electron Green function in magnetic field 
has been investigated by Ando and Uemura |]2T[] and in more detail by Laikhtman and E. 



Altshuler ||22[ for Fg determined by an elastic short range scattering. In the case we are 
particularly interested in, i.e. GaAs the temperature variation of mobility from the liquid 
helium temperature to the temperature of experiment (about 200K) is substantial (in the 
typical cases, at least, by several times). It means that the acoustic phonon scattering 
(that can be considered as short range elastic) should be predominant. It was shown in 



Refs. PT,22| that in this case the electron Green function has a non-Lorentzian form with 



the characteristic width Fg given by 



Fe = F(°)vW^, F(°) = ^n/2TTT (3.4) 

where r is the relaxation time for B = obtained by assuming the same scatterers as for a 
finite B. 

Further on our formulas should be considered as order-of-magnitude estimates giving the 
correct dependencies on the parameters, though not the parameter-independent numerical 

11 



coefficients of the order of unity. For such estimates it will be sufficient to use the Lorentzian 
form of n(u;,q). As has been indicated, in the resonance approximation one should retain 
only the resonant term of all the series ( p.lO| ) for Il{i{uj, q) 



Ur{uj, q = -r- 3.5 

UJ — TV cosy + tl e 

where TZj^ is the residue at the pole u = ATi^cos^^ — iT^. Calculating a^x one can evaluate 
the integral over the frequencies taking the residues in the poles uj = N'QcosO ± iFg. To 
get the result one should remove the factor S{6~^) in Eq. ( |2.20|) and insert instead into the 
integrand of Eq. ( p.21| ) 

Img;^^ + 27 



where 



Fe 27re^ 

7 = ^; LJ = 7l_\f{AfQcos6,q) 

UJ q 



while Ea should be calculated at = N'Q cos + 2(F + Fg). One can see that integral ( |2.21| ) 
is dominated by the values of q where the asymptotic expansion of the Bessel function is 
valid, so that 



2 



We remind the reader that Eq. (|2.2CI|) is derived within the so-called RPA (loop) ap- 



proximation [see Eq. ( p.l6|) ] as it takes into account the resonant interaction of electrons 



with optic phonons. Such resonant terms should probably be added also to the vertex parts 
describing the polaron effects. These terms may be essential in the magnetophonon maxima. 
However, we are looking for the critical values of the parameters such as electron concentra- 
tion where the resonant interaction disappears, so that the MPR signal sharply goes down. 
In such a situation the polaron effects should be also suppressed. Thus we believe that our 
theory, though probably not permitting the exact calculation of the MPR amplitude, still 
gives the correct characteristic values of the electron concentrations and the angles where 
the MPR signal rapidly decreases. 
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In the present paper we limit ourselves with the sufficiently low concentrations Ug where 
the MPR maxima are well defined. As is shown in Ref. [|10|, an additional mechanism of 
electron level broadening due to the electron-electron interaction appears at high electron 



concentration [see Eq. (6.10) of Ref. [|T0|)]- Here we consider the concentrations that are 
not so large that one would have to take into consideration this effect. 

For small values of 7 Eq. (|3.6| ) turns into \m.eR{uJt) as one can neglect the terms 27. 
When 7 has reached the critical value (l/2)|£^^(c<Jt)| the amplitude of the maximum begins 
to go down. This condition can be written as 

27 ^ Im {eaY^ . (3.8) 

for 

NVLcose = ujf (3.9) 

Under the MPR condition ( p.9| ) Im {£a)~^ is small. Therefore even a small variation of the 
quantities determining 7 can violate Eq. (p.8|). The screening is important provided that 



the parameter 

(3 = lm [e a{MVL cos e + iVe + iT)]~^ (3.10) 

is bigger than (or of the order of) unity. If it is much bigger than unity the screening is 
strong and the result is independent of /3. The MPR maxima become /3-dependent when (5 
is of the order of unity. Then a relatively small variation of (3 might drastically change the 
result (see Figs. 4 and 5 of Ref. [0). To achieve such a change one needs a variation of 
that need not be large. 

Let us follow a variation of the MPR signal as a function of 9. One can write for an 
arbitrary magnetic field B{9) that is near the A/"th resonance value Bj^{9) 

B{9) = B^{9) + AB. (3.11) 

-^^"(0) = huJtmc/eAf is the position of the A/'th MPR maximum for ^ = 0. The field 
corresponding to the A/'th resonance, varies for finite values of 9 according to 
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Bm{9) COS 9 = B^{0). (3.12) 

Thus the MPR amphtudes for various angles 9 (and the same resonance number JV) corre- 
spond to the same value of jS and therefore should coincide. (Only the width of the maximum 
should enhance with decrease of 9) . This conclusion of the theory is in drastic disagreement 
with the experiment 0. 

We think that this is due to the assumption that the electron concentration is a 
constant independent of B. It is known, however, that in the course of temperature variation 
from 77 K to 300 K at i? = the variation of the carrier concentration may comprise several 
tens per cent. This means that the energy variation of some donors on the scale of the order 
of hundred K, or so noticeably shifts the electron concentration balance between the donors 
and the well. 

Let us give a rough estimate of the variation of impurity level positions as a function of 
B. If one assumes a hydrogen-like states the variation 6e of their lowest level position e{B) 



under the shift of magnetic field 6B is (for instance, [|T4|, section 112) 



B 2\n[hn{B)/eo] ^ ^ 

where eo = me^/2/i^e^ is the Bohr energy while e{B) is the position of the level in magnetic 
field B, mark that e{B) < 0. This is a lower estimate as with the overlap of atomic orbits 
the influence of the magnetic field should enhance. As in our case ln[hQ{B) / e^] is, roughly, 
of the order of unity we will accept the estimate 6e/e ~ —6B/B. This should result in 
redistribution of electrons between the donors and the well. In other words, the electron 
concentration in the well will decrease (see below). The chemical potential /x can be obtained 
from the equation 



Ni = exp ( ) { Ni exp 



S 

2^ exp 



27ra^ cos^ 9 ^ 



hQiiB cos 9) 



(3.14) 



where Ni is the total number of donors. Further on we will assume that the number of 
electrons bound to the donors is much bigger than in the well. This seems to be a rather 
typical experimental situation. Practically the actual form of the impurity states and the 
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distribution of impurity levels may be (and usually is) much more complicated. What is 
actually relevant to bring about a variation of the electron concentration Ug in the well is 
the magnetic field dependence of the positions of impurity levels in the proper interval of 
energies. 

The variation of electron concentration for the values of B given by Eq. ( |3.12| ) is 

with ^ (3.15) 



ns{B) ^B l + ns{B)/ni{B) 

where ni{B) is the 2D concentration of the electrons bound by the donors. It will be natural 
to assume below for the estimates that |.^| is of the order of unity. 

It is convenient to present /3 as a ratio of the electron concentrations to some characteristic 
value 

P = ns{B)/ndown (3.16) 

where 

1 V2^he^uJtB^{0) Te + r ^g^^^ 



^down SoaQxT'^^emC UJl - UJt 

^down is the lower critical concentration where the screening ceases to play a role. For 7V=3 
in GaAs its characteristic value is about 10^° cm~^. Equation (|3.17|) differs from the equation 



for ridown given in Ref. [|T0| as we have assumed there qt ~ and hut ~ 2A;bT that is valid 



for GaAs under certain conditions and may differ for other situations. 

The parameter (3 depends on the electron concentration n^. This fact permits one to 
compare the dependence of the MPR amplitudes in perpendicular magnetic field under the 
variation of concentration (i.e., in different samples) and the angular dependence of the 
MPR maxima under variation of 6. This will permit to check as to whether the magnetic 
field induced variation of is sufficient to explain the decrease of the MPR amplitude as 
a function of 6 on the one hand and the existence of the subsidiary (side) maximum in the 
angular dependence of the MPR amplitude on the other hand (see Ref. 0). 

We start with discussion of the behavior of the MPR amplitude under magnetic field 
rotation at small angles 0^1. We will consider the concentration interval Ug ^ rido 



aown 
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where there is a sharp dependence of the MPR amphtude on p. According to Eq. ( |3.12|) the 
MPR maximum will shift by 

5Bm^ B^{0){1~ cos 9). (3.18) 

This will result in the relative variation of concentration Sus/ng ^ ^(1 — cos 6) that manifests 
itself in the angular dependence of the MPR amplitude. As is known from the experiment 
in the perpendicular magnetic field, in this region of concentrations when goes down it 
brings about a sharp decease of the amplitude. 

As a reasoning supporting our view we will consider the following numerical example. 
According to Ref . P , Fig. 6 the variation of the MPR amplitude is about 60 % provided 
the concentration varies within the interval from 1.8-10^'' cm~^ to 2-10^'' cm~^. The same 
decrease of the amplitude induced by magnetic field should take place for 5ns/ ~ 0.1, i.e. 
for 9 of the order of 20°. This shows an order-of-magnitude correspondence with the results 
of Ref. . Further increase of 9 gets the system into the region where n < ridown and a well 
defined resonant peak with the resonant condition (|3.9| ) disappears. 

The behavior of this sort takes place provided that 1. For such concentrations where 
/3 ^ 1 small variation of this parameter does not play an essential role. This means that 
for the samples that have a maximal MPR amplitude (n^ ~ 5-10^° cm~^) its 9 dependence 
should be absent. 

Now we will treat the region of large angles. Consider the samples with the concentration 
of the order of ndown- The angle 9 going up, can become so small that the perturbation 
approach becomes applicable. In other words, one observes the resonance determined by 
the condition 

MVLcos9 = uJi. (3.19) 



The screening is not important provided that (3 <1. In view of Eq. (p.l9|) one should insert 
in (|3.10|) SAioJi + iVe + iV). As a result, one can write 



A = ^-f^V. (3.20) 



^down V Fg ~l~ F 
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If the last factor in this equation is large enough one needs rather small concentrations ns{6). 
Let us estimate the angles where they can be achieved. It follows from Eq. ( |3.14| ) under the 
same assumptions provided that 6B/B is not small 

nsi6B + B)= UsiB) exp{6e/kBT). (3.21) 



Very rough estimates in the spirit of Ref. ||26| give [the dependence e{B) is rather smooth 



and we believe that its expansion up to the linear term can be justified] 

Se ^ -Ceo (^^^ = -C^^ (3.22) 

where C is a number of the order of unity. This estimate is based on the idea that the level 
shift depends on a single parameter, i.e. the ratio of the magnetic energy to the Coulomb 



one (cf. the analytical treatment of the hydrogen atom in magnetic field in Ref. [|1^, §112). 



Thus the discussed cj^-niaximum exists for the angles bigger than 6c given by the equation 



The second term under the sign of root is, most probably, of the order of unity. In this case 
6 is somewhere in the interval between 40° and 70°, or so. Thus in the regions of small and 
large angles we have well-defined ut and cuz-resonances, respectively. Their positions depend 
on BcosO. In the intermediate interval of angles there are no well-defined MPR maxima. 

It is interesting to note that in the region of small angles the width of a maximum goes 
up with the number of resonance A/". It can be easily seen if one takes into consideration 
enhancement of with decrease of magnetic field as well as enhancement of ric with A/". All 
these conclusions are in a qualitative agreement with the results described in Ref. p. 

IV. CONCLUSION 

To summarize, we would like to stress that the interpretation of behavior of the MPR 
in a tilted magnetic field has been a long-standing problem Two types of resonant 
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maxima have been discovered on experiment, i.e. the oot- and turresonances as their positions 
are determined by the frequencies Ut and ui, respectively. 

Important points to provide theoretical interpretation of these resonances are the de- 
pendence of all quantities characterizing the 2D motion of electrons on the combination 
BcosO whereas their total concentration in the well depends on B. Had only the BcosO 
dependence existed, the amplitude of the resonances, i.e. the heights of the maxima would 
have been independent of the angle 6. We think that the only way to preserve the theory of 
2D electron gas describing the concentration dependence of the MPR is to assume variation 
of the electron concentration with B. It manifests itself in the angular dependence of the 
MPR amplitude. In this way we have been able to give a qualitative interpretation of the 
results given in Ref. 0. 
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APPENDIX: CALCULATION OF POLARIZATION OPERATOR 

As indicated in Sec. ||, we assume for the electrons a parabolic confining potential 
mujQz'^/2 with the gauge for the vector potential A = {—By cos6 + 5^ sin 6', 0, 0). In spite 
of the first inequality ( p.8| ), it is convenient to look for the exact transformation of the 
Hamiltonian and solution of the Schrodinger equation and only then go to the limit 



Applying a standard procedure of diagonalization we get for the bigger eigenvalue 



n/uo < 1. 



(Al) 



(A2) 



while the smaller eigenvalue is 
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2 _ + ^0 



1 - 



(A3) 



The variables (Y, Z) diagonalizing the Hamiltonian are expressible through the initial 
variables (j/, z) as 



Here C = 1/VTT 



}_( y-yo-rz \ 

C \r{y -yo) + zJ ' 

m - cos2 e 



(A4) 



^2 COS 6 sin ^ 

As a result, we get two noninteracting oscillators, i.e. the F-oscillator and the Z-oscillator 

1 92 1 ^2 



2m dY^ 2m dZ^ 2 



(A5) 



We are interested in the eigenfrequency of the y-oscillator that is 



= ^ cos 9 



(A6) 



plus small terms proportional to Vl^ /uJq <^ 1. The eigenfrequency of the Z-oscillator is equal 
to ojq (with the same accuracy). 

Now we turn to calculation of the polarization operator for nondegenerate free elec- 
trons in magnetic field by a slight modification of the method proposed by Sondheimer and 



Wilson [^. The method is based on the spectral representation (see Sec. 36 [0). 



^ / / N {m\niQ , z')\l) {l\niQ , z)\m) ,^ , , . „, 

T{r{u,z',z) = Y^w r ' ^ ' ^' ' V ■ [1 -exp(-^i„/T)]. (A7) 

Sondheimer and Wilson introduced a Green function of a complex time argument 7 



G(r, r', 7) = E exp {-Sp-f) , Re 7 > 0. 



(A8) 



Here (3 is the set of all quantum numbers of an electron, while ep is the energy of the 
quantum state (3. One can express the polarization operator through the Green functions 
of a complex time argument. 

The retarded polarization operator for Hamiltonian ([A5|) is 
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ng)(r, r', t) = 2^ cosh exp (^) e(t) 



X 



G ( r', r, i - zt) G(r, r', zt) - G{r', r, (r, r', ^ + it 



(A9) 



Here is a product of the eigenf unctions of the Y- and Z-oscillators, 0(t) is the step 

function. 

In order to sum up the series in Eq. over the Y- and Z-oscillator quantum numbers 



we will use the following relation p 



OO 



1 - e 



-1/2 



exp 



A2l + e-^" 



1 1 — e""^*^ 1 



(AlO) 



where Hn{y) is the Hermite polynomial. It is convenient to replace the summation over the 
quantum number by integration over y^. 
Now, 

Gfr, r', 7) = "^Q^ f sinh 'jilf, sinh 7i7n^ f tanh — — - + tanh — — ^ ^ 



xexp{4coth(2^)(Ar)^-i|coth(2^)(AZ)^-l5?tanh(^)( 
1 C^lAx + ib tanh {no-f/ 2) {z + z')f \ 



z + z 



i\2 



4 tanh {-iVtj2) + {a{r^/ai) tanh (7ilo/2) 



(All) 



Here 



.air b 
o = — , oi = -; 
^2 r 



the variables in Eq. (|A11|) are made dimensionless by the transformations 



(A12) 



AY = (Y -Y')/au AZ = {Z - Z')/ai, z ^ z/ai, z'-^z'/a^ 



(A13) 



where 



(AM) 
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To calculate the polarization operator it is convenient to use the Green functions in the 
momentum representation. After rather involved but straightforward calculations we get 



n^'Hq, 9., Qz. t) = p(q, g„ g„ t) + p*(-q, g„ t) 



(A15) 



where the asterisk means a complex conjugation; 



'2 , ^^(ll 

Hz ' 



2 



2iCaiB'^ cosh 



2T 



7ra2 sinh(^],/2^) sinh(^]o/2T) ' ""^^ [ + (5 ) ^^^^^^ 



C2 



af 2 



(AI6) 



where 



^?2(^^) = I 



coth — - — cos \\lt -\ — — 
2T V 2T 

O in 

coth — - + cos [\lt -\ — — 
2T V 2T 





■ 


sinh 




2T. 




1] ■ 


sinh 






2T. 



sin^](t + i/2^) 



sinh(fi/2T) 

Here q and refer to the Fourier components over the differences r — r' and z — z' while 
Qz is related to {z + z')/2. 

We are looking for a frequency representation of the polarization operator. We are going 
to take into consideration only the lowest level of transverse quantization, i.e. the lowest 
miniband. At the first sight one could try to average Eq. ( |A16| ) over the time interval bigger 
than 1/ujq but smaller than l/ficos^^. However, some spurious terms can appear as a result 
of the direct averaging due to the even powers of sin uqI and cos Uot. Therefore we will look 
for the frequency representation of the whole expression Eq. (|A16|) . One can find it with 
the help of the identities 



exp(-2:cosy4) = ^e™^/„(z), exp(-2: sin A) = ^ e*"^ J„ 



iz] 



(A17) 



where Jm are the ordinary Bessel functions. Combining them with Eq. ( [A16| ) one can see 
that y4 is a linear function of time, so that one can easily calculate the Fourier components. 
Discarding the poles describing the transitions to the higher minibands one gets 
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= -2Cns exp 



2 2C2 



X 



2T C2 \^8af 
"+'1^2sinh(fi,/2T)y '-"1^8C2sinh(fi,/2T) 



+ — coth 



2T 



a? I 2 ^ 2C2 i 



Imn 
X J, 



irqxQz 



sinh(ficn/2T) 



(A18) 



2^2 sinh(^],/2^) / - fi^^ + z5 ' 

There are elegant formulas to calculate this sum. However, as we are interested in the 
limiting case Q'^/uq <^ 1 one can see the result immediately from Eq. (|A18|) . Indeed, the 
Bessel functions having in their arguments as a factor the small parameter r can be discarded 
unless I — m = m = and we are left with 



n 



(3) 
R 



= —2ns exp 



2 ^° 2T al 8 a\2 



+ 



sinh {nnj2T) 



(A19) 



2 sinh(fic/2T) y uj — Vtji + id 

One can see that H^^) Eq. (|MUD coincides with ng^ of Eq. and if one takes 

into consideration that within the accepted approximation VLc = cos and that we use in 
Eq. ( [A19| ) dimensionless variables. 
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